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Abstract 

We investigate the status of the neutrino decay solution to the solar neutrino problem 
in the context of the recent results from Sudbury Neutrino Observatory (SNO). We present 
the results of global x^-analysis for both two and three generation cases with one of the 
mass states being allowed to decay and include the effect of both decay and mixing. We find 
that the Large Mixing Angle (LMA) region which is the currently favoured solution of the 
solar neutrino problem is affected significantly by decay. We present the allowed areas in the 
Am^ — tan^ 6 plane for different allowed values of a and examine how these areas change 
with the inclusion of decay. We obtain bounds on the decay constant a in this region which 
implies a rest frame life time To/m2 > 8.7 x 10~^ sec/eV for the unstable neutrino state. 
We conclude that the arrival of the neutral current results from SNO further disfavors the 
neutrino decay solution to the solar neutrino problem leaving a very small window for the 
decay constant a which could still be allowed. 
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1 Introduction 



The comparison of the SNO charged current (CC) measurments with their neutral current (NC) 
rate estabhshes oscillations to active neutrino flavor at the 5.3a level [|l|. The inclusion of Super- 
Kamiokande (SK) electron scattering (ES) rate with SNO confirms transitions to active flavors 
at more than 5.5a level. In other words SNO rules out the possibility of having a purely sterile 
component in the solar neutrino beam at the 5.5a level. However presence of sterile components 
in the beam is not entirely ruled out. Transition to mixed statesQ is still allowed with < 32%(65%) 
sterile mixture at the la{2a) level H, H. One of the scenarios in which one can have a final 
sterile state is neutrino decay. The aim of this paper is to investigate the status of the neutrino 
decay solution to the solar neutrino problem with the incorporation of the SNO data. We deal 
with both two generation and three generation cases. For the latter we consider three flavor mixing 
between z/g, z/^ and u^- with mass eigenstates i>i, V2 and z/3 with the assumed mass hierarchy as 
m\ < 1712 < "^3- The lightest neutrino mass state is assumed to have lifetime much greater than 
the Sun-Earth transit time and hence can be taken as stable. Whereas the heavier mass states z/2 
and z/3 may be unstable. But the CHOOZ data |^ constrains the mixing matrix element Ues to a 
very small value which implies very small mixture of the z/3 state to the z/g state produced in Sun. 
So the instability of the z/3 state has hardly any effect on solar neutrino survival probability, and 
to study the effect of decay on solar neutrinos for simplicity we can consider only the second mass 
state z/2 to be unstable. 

There are two kinds of models of non-radiative decays of neutrinos 

• If neutrinos are Dirac particles one has the decay channel 1^2 — > z/i/j -|- (p, where is a right 
handed singlet and is an iso-singlet scaler |^. Thus all the final state particles are sterile. 

• If neutrinos are Majorana particles, the decay mode is z/2 ^ z/i + J, where z/i interacts as a 
z/e with a probability |?7eiP and J is a Majoron |Q. 

In both scenarios the rest frame lifetime of z/2 is given by 



where g is the coupling constant and Am^(= — mf) is the mass squared difference between the 
states involved in the decay process. If a neutrino of energy E decays while traversing a distance 
L then the decay term exp{-aL/E) gives the fraction of neutrinos that decay, a is the decay 
constant and is related to tq as a = ^2/70. 

For the Sun-Earth distance of 1.5 x 10^^ m, and for a typical neutrino energy of 10 MeV one 
starts getting appreciable decay for a ~ 10~^^ eV^. For lower values of a the exp{-aL/ E) term 
goes to 1 signifying no decay while for a > 10~^° eV^ the exponential term goes to zero signifying 
complete decay of the unstable neutrinos. Assuming m2 » mi the equation (0) can be written 



IQtt 7712(1 + 1711/7712) 
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^0 = 



(1) 



as 



(2) 



^States which are a mixture of active as well as sterile components. 
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If we now incorporate the bound < 4.5 x 10 ^ as obtained from K decay modes |]T0| we get the 
bound Am^ > lO^a. This gives the range Am^ ~ 10~^ eV^ for which we can have appreciable 
decay. 

Neutrino decay in the context of the solar neutrino problem has been considered earher in 
[0, ^, |11[] and more recently in |]12|, |T^, |1^]. After the declaration of the SNO CC data last year 
p!5[l , bounds on the decay constant a have been obtained from the data on total rates [0 as 
well as from SK spectrum data [T^. However with the declaration of the recent SNO neutral 
current results [0] we expect to have a better handle on the sterile admixture in the solar neutrino 
beam and hence tighter constraints on the decay constant a. The statistical analysis for the two- 
generation oscillation solution for stable neutrinos with the global solar neutrino data including 
the recent results of the SNO experiment can be found in f^, |18], |19|. In this paper we consider the 
possibility for unstable neutrino states and do a combined global satistical analysis of the solar 
neutrino data including the total rates from the CI and Ga experiments pO|, the zenith-angle 
recoil energy spectrum from SK [|2[ and the combined (CC+ES+NC) day-night energy spectrum 
from SNO [||]. We do our analysis for both two and three generation scenarios. We include the 
results of the CHOOZ experiment [^] in our three generation analysis. 

In section 2 we briefly present the survival probability for the solar neutrinos with one of 
the components unstable. In section 3 we first find the constraints on a coming from the global 
solar neutrino data within a two-generation framework. We argue that decay of solar neutrinos 
is largely disfavored as a result of the inclusion of the SNO day-night spectrum. We next extend 
our analysis to include the third neutrino flavor and present bounds on a, Am^i and tan^ 621 for 
different allowed values of ^13, from a analysis which includes both the global solar neutrino 
data as well as the CHOOZ data. We end with conclusions in section 4. 



2 Formalism 

The three-generation mixing matrix relating the mass and flavour eigenstates are given as 
U = -R23-R13-R12 

^ C13C12 ■S12C13 Si3 \ 

= — S12C23 — ■523'5l3Ci2 C23C12 — S23S13S12 S23C13 (3) 
\ ■523S12 — S13C23C12 — S23C12 — ■513S12C23 C23C13 / 

where we neglect the CP violation phases. 

Allowing for the possibility of decay of the second mass eigenstate the probability of getting a 
neutrino of flavour / starting from an initial electron neutrino flavour is 0. 

-Te/ — Ctel l^l/l + '^e2 l^2/| 6 + «e3 1^3/1 

where Ei is the energy of the mass eigenstate i, E is the energy of the neutrino beam, a is the 
decay constant, L is the distance from the center of the Sun, Rq is the radius of the Sun, 0j_0 are 

^For a rigorous derivation of the probability for the decay plus oscillation scenario see 12^. 
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the phase inside the Sun and a% is the amphtude of an electron state to be in the mass eigenstate 



Ui at the surface of the Sun 22 



a- = y: x,,uf: (5) 



,0^ _ 

i=l,2,3 



where Xij denotes the non-adiabatic jump probabihty between the i and j states inside the Sun 
and f/j® denotes the mixing matrix element between the flavour state and the mass eigenstate 
Pj in Sun. Afr denote the z/j — > Uf transition amplitudes inside the Earth. We evaluate these 



amplitudes numerically by assuming the Earth to consist of two constant density slabs |2^. It 
can be shown that the square bracketed term containing the phases averages out to zero in the 
range of Am^ in which we are interested [^. For further details of the calculation of the survival 



probability we refer to p2 



In the two generation limit the Ve survival probability including Earth matter effects and decay 



can be expressed as [|13 



p _pday I (sin^g-P,,)(P,t^ + e-^"(^-^e)/^(P,t^-l)) 

where is the transition probability of z/2 — ^ v^. at the detector while P^^^ is the day-time 
survival probability (without Earth matter effects) for v^. and is given by 



Pf:^ = cosH[Pjsm^eM + il-Pj)cos^OM] 

+ sin^ ^[(1 - Pj) sin^ Om + cos^ ^M^jje^'^^^"^®)/^ (7) 

where Pj is the non-adiabatic level jumping probability between the two mass eigenstates for 
which we use the standard expression from |23] and 9m is the matter mixing angle given by 



Am^ sin 29 , , 

tan2^A/ = ^ • 8) 

Am2 cos 20 - 2y2Gi.neE 

Tie being the ambient electron density, E the neutrino energy, and Am^ (= m| — m^) the mass 
squared difference in vacuum. 

From eq. (0) we note that the decay term appears with a sin^ 9 and is therefore appreciable 
only for large enough 9. Thus we expect the effect of decay to be maximum in the LMA region. 
This can also be understood as follows. The i^e are produced mostly as V2 in the solar core. In the 
LMA region the neutrinos move adiabatically through the Sun and emerge as V2 which eventually 
decays. For the SMA region on the other hand Pj is non-zero and Vf. produced as z/2 cross over 
to z/i at the resonance and come out as a from the solar surface. Since is stable, decay does 
not affect this region. 



3 Analysis of data and results 

3.1 Bounds from two-generation analysis 

First we use the standard minimisation procedure and determine the Xmm and the best-fit 
values of the oscillation and decay parameters, in a two flavor mixing scenario. Details of our 
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statistical analysis procedure, including the definition of the and the correlated error matrix 
are given in Appendix A. We incorporate the total rate in CI and Ga, the zenith-angle spectrum 
data of SK and the SNO day-night spectrum. We find that the global best-fit| comes in the LMA 
region with the decay constant a = and tan^ 6 = 0.41, Am^ = 6.06 x 10~^eV^ . 

In figure 1 we plot the Ax^ (= ~ xLin) ^^e two generation neutrino decay scenario, 
against the decay constant a, keeping Am^ and tan^ 6 free. The figure shows that the fit becomes 
worse with increasing value of the decay constant. We remind the reader that neutrino decay 
is important only for LMA and for each value of a on this curve the minimum comes with 
Am^ and tan^ 6 in the LMA region. The inclusion of the SNO data in the global solar neutrino 
analysis improves the fit in the LMA region for stable neutrinos^. Hence for small values of a 
where neutrinos can be taken as almost stable, the inclusion of the SNO data gives a better fit. 
However as a increases the neutrino decay leads to more sterile components in the final state - 
which is disfavoured by the SNO/SK combination and the fit worsens with SNO included. As is 
seen from figure 1 the global solar neutrino data put an upper bound on the decay constant - 
a < 7.55 X 10~^^eV^ at 99% C.L., when Am^ and 6 are allowed to take on any value. Before the 
declaration of the SNO results the bound on a from combined analysis of CI, Ga and SK data 
was a < 3.5 x 10~^^ eV^ |T^. Thus the inclusion of the recent SNO results have further tightened 
the noose on the fraction of neutrinos decaying on their way from the Sun to Earth and neutrino 
decay is now barely allowed. 

In figure 2 we plot the allowed areas in the Am^ — tan^ 6 plane for different allowed values of the 
decay constant a for a two flavor scenario and with global data including SNO. The figure shows 
that the allowed area in the LMA region is reduced as a increases. As a increases the spectral 
distortion increases and the effect is more for higher(lower) values of tan^ ^(Am^) |]13|. Since both 
the SK as well as the SNO spectra are consistent with no energy distortion, these regions get 
disallowed with increasing a. Spectral distortion due to decay is less for the high Am^. However 
these values for the Am^ are now disfavored at more than 3a by SNO since the values for the 
flux required to explain the suppression in SK and SNO CC measurements are much lower than 
those consistent with the SNO NC observation. In addition the low energy neutrinos relevant for 
Ga experiment decay more making the fit to the total rates worse. 

As discussed earlier, for small values of mixing angles the fraction of the decaying component 
in the z/g beam (~ sin^ 9) being very small decay does not have much effect. Hence no new feature 
is introduced in the SMA region because of decay and it remains disallowed by the global data. 
On the other hand, in the LOW region though mixing is large, due to small Am^ any appreciable 
decay over the Sun-Earth distance is not obtained, if the bound on the coupling constant from 
K-decay is to be accounted for. Therefore, LOW region in the parameter space remains 
unaffected by the unstable z/2 state, and we do not show it explicitly in figure 2. 



■^The bounds that we give here and subsequently apply to the decay model where both final states are sterile. 
In the Majoron decay model the final state 17^ can interact in the SK/SNO detectors. However there is an energy 
degradation and the best-fit values do not change significantly An interesting possibility where the absolute 
mass scales of the two neutrino states are approximately degenerate and hence the daughter De is not degraded in 
energy is considered recently in . 

^See [||, |l9|, ^ for recent global analyses and |26, 2^ for earlier ones. 
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3.2 Bounds from three-generation analysis 



Next we investigate the impact of the unstable 1^2 mass state on the three flavor oscillation scenario 
with Am2i = AttIq and Arrig-^ = Am^^i^QQ^ ~ Am^^^ = Amgg. For the three generation case the 
is defined as 

(9) 

The x^soiar is defined in Appendix A while for the expression of Xc/1002 '^^ refer the reader to 
[p^ . In figure 3 we plot the Ax^(= — Xmin) ^ various fixed values of 61^ allowing the 
other parameters to take arbitrary values. Am|]^ is varied within the range [1.5 - 6.0] x 10"'^ eV^ 
as obtained from a combined analysis of atmospheric and CHOOZ data This figure 

indicates the allowed range of a for a fixed 6*13. As ^13 increases the contribution from xl.hooz 
increases shifting the plots upwards. Consequently the curve corresponding to higher tan^ 6'i3 
crosses the 99% CL. limit at a lower a and a more stringent bound is obtained as compared to 
smaller tan^ cases. 

In figure 4 we plot the allowed areas in the tan^ 9i2 — Am^i plane for different sets of values of 
tan^ ^13 and Amg^^ at different values of a. For fixed values of tan^ ^13 and Amgj^ the allowed area 
decreases with increasing a because of increased spectral distortion and transition of the solar 
neutrinos to sterile components. At any given a the allowed area shrinks with the increase of 
tan^ 6*13 and Amg^. This is again an effect of the CHOOZ data. 



4 Conclusions 

In this paper we have done a global x^ analysis of the solar neutrino data assuming neutrino 
decay. We incorporate the full SNO (CC+ES+NC) day-night spectrum data. We present results 
for both two and three generation scenarios. We find that the best fit is obtained in the LMA 



region with the decay constant a zero i.e. for the no decay case. We had pointed earlier |T2|, |13 
that the decay scenario was in conflict with the data on total rates because it suppresses low energy 
neutrinos more than the high energy ones. However decay also predicts distortion of the neutrino 
energy spectrum and presence of sterile components in the resultant neutrino beam, both of which 
are severely disfavored by the recent results from SNO and SK. For the two generation case the 
bound at 99% C.L. on a that one gets after including the SNO data is a < 7.55 x 10~^^eV^. This 
corresponds to a rest frame lifetime tqIvti^ > 8.7 x 10~^s/eV. This is consistent with astrophysics 
and cosmology. The effect of decay is not important in the SMA region because the decaying 
fraction in the solar z/g beam goes as sin^ 9 and this region remains disallowed. There is no change 
in the allowed areas in the LOW region also because of eq. (0) and the bounds on the coupling 
constant from K decays which restricts Am^ ~ 10~^ eV^ and decay over the Sun-Earth distance 
does not take place in the LOW region. The three generation analysis gives stronger bounds on 
a for non-zero ^13. In conclusion, the inclusion of the recent SNO data severely constrains decay 
of the solar neutrinos such that they are largely disfavored. However nonzero values of a remain 
allowed at the 99% C.L.. 



^The range of '^m'^i allowed at 99%C.L. from the combined analysis of final SK+K2K data is [2.3 - 3.1] x 10 
eV^ [pOj. However the values of Amg^^ that we use are still allowed at the 3a level. 
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A The statistical analysis of the global solar data 



We use the data on total rate from the CI experiment, the combined rate from the Ga experiments 
(SAGE+GALLEX+GNO), the 1496 day data on the SK zenith angle energy spectrum and the 
combined SNO day-night spectrum data. We define the function in the "covariance" approach 
as 

N 

where is the number of data points (2 + 44 + 34 = 80 in our case) and {cfj)''^ is the inverse 
of the covariance matrix, containing the squares of the correlated and uncorrelated experimental 
and theoretical errors. The only correlated error between the total rates of CI and Ga, the SK 
zenith-energy spectrum and the SNO day-night spectrum data is the theoretical uncertainty in 
the flux. However we choose to keep the flux normalization Jb a free parameter in the 
theory, to be fixed by the neutral current contribution to the SNO spectrum. We can then block 
diagonalise the covariance matrix and write the ^is a sum of for the rates, the SK spectrum 
and SNO spectrum. 

2—2,2,2 (T'W 
X — Xrates + XsKspec + XsNOspec \^^) 

For xLes we use i?cT = 2.56±0.23 SNU and ii^f = 70.8±4.4 SNU. The details of the theoretical 



errors and their correlations that we use can be found in |2^, |3T . 

For the 44 bin SK zenith angle energy spectra we use the data and experimental errors given 
in 0. SK divides it's systematic errors into "uncorrelated" and "correlated" systematic errors. 
We take the "uncorrelated" systematic errors to be uncorrelated in energy but fully correlated in 
zenith angle. The "correlated" systematic errors, which are fully correlated in energy and zenith 
angle, include the error in the ^B spectrum shape, the error in the energy resolution function and 
the error in the absolute energy scale. For each set of theoretical values for Am^, tan^ 6 and a we 
evaluate these systematic errors taking into account the relative signs between the different errors 
[ |T8|1 . Finally we take an overall extra systematic error of 2.75% fully correlated in all the bins [§. 

Since it is not yet possible to identify the ES, CC and NC events separately in SNO, the SNO 
collaboration have made available their results as a combined CC+ES-I-NC data in 17 day and 17 
night energy bins. For the null oscillation case they do give the CC and NC (and ES) rates 
but firstly these rates would slightly change with the distortion of the ^B neutrino spectrum from 
the Sun and secondly they are highly correlated. They work very well for studying theories with 
little or no energy distortion such as the LMA and LOW MSW solutions if the (anti) correlations 
between the CC and NC rates are taken into account and can be used to study the impact of 
the NC rate of SNO on the oscillation parameter space ||^. They can also be used to give insight 
into the ^B flux and it's suppression by analysing the SK and SNO CC and NC data in a model 
independent way 0. However for theories like neutrino decay where we expect large spectral 
distortions they cannot be used. Hence we analyse the full day-night spectrum data in this paper 
by adding contributions from CC, ES and NC and comparing with the experimental results given 
in [jl| . For the correlated spectrum errors and the construction of the covariance matrix we follow 



the method of "forward fitting" of the SNO collaboration detailed in [32 
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Figure 1: Ax {= X ~ Xmin) decay constant a, for two-generations. Also shown is the 99% 
CL. limit for three-parameters. 
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Figure 2: The 90, 95, 99 and 99.73% C.L. allowed area from the two-generation global analysis of 
the total rates in CI and Ga, the SK zenith angle energy spectrum data and the SNO day-night 
spectrum in presence of decay and oscillation. The C.L. contours are for Ax^ corresponding to 
three parameters. 
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Figure 3: The Ax^(= — Xmin) decay constant a for the three generation case for various 
fixed values of 6*13. The other parameters are allowed to take any value. 
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Figure 4: The 90, 95, 99 and 99.73% CL. allowed area from the three-generation analysis of the 
global data in presence of decay and oscillation. 
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